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A GENERALIZATION OF THE DIGITAL BINOMIAL THEOREM 


HIEU D. NGUYEN 


Abstract. We prove a generalization of the digital binomial theorem by constructing a one-parameter 
subgroup of generalized Sierpinski matrices. In addition, we derive new formulas for the coefficients of 
Prouhet-Thue-Morse polynomials and describe group relations satisfied by generating matrices defined in 
terms of these Sierpinski matrices. 


1. Introduction 


The classical binomial theorem describes the expansion of {x + y)^ in terms of binomial coefficients, 
namely for any non-negative integer we have 

N / \ 

( 1 ) 

where (^) are defined in terms of factorials: 

/A\ m 

\k) ~ k\{N-k)\' 

In [^, the author introduced a digital version of this theorem where the exponents appearing in m are 
viewed as sums of digits. To illustrate this, consider the binomial theorem for N = 2: 

(x + y)^ = x^y° + x^y^ + x^y^ + x^y^. ( 2 ) 

It is easy to verify that m is equivalent to 

{x + = ^Si3)ysi0) ^ ^si 2 )ySil) ^ ^S(l)y.i 2 ) ^ ^s(0)ySi3)^ (3) 

where s{k) denotes the sum of digits of k expressed in binary. For example, s(3) = s(I • 2^ + I • 2°) = 2. 
More generally, we have 


Theorem 1 (Digital Binomial Theorem [^). Let n he a non-negative integer. Then 

(a: + 2/)®(”)= ^ j.s{rr^)ys(n-m) 


0<m<n 

{m,n—m) cari^-free 


(4) 


Here, a pair of non-negative integers (j, k) is said to be carry-free if their addition involves no carries when 
performed in binary. For example, (8, 2) is carry-free since 8 -I- 2 = (1 • 2^ -1- 0 • 2^ -|- 0 • 2^ -|- 0 • 2°) -|- I • 2^ = 10 
involves no carries. It is clear that (j, k) is carry-free if and only if s(j) -|- s(fc) = s(j -|- k) (see [2112]. Also, 
observe that if n = 2^ — 1, then (Hj) reduces to (|T]). 

In this paper we generalize Theorem [I] to any base b >2. To state this result, we shall need to introduce 
a digital dominance partial order on N as defined by Ball, Edgar, and Juda in [2]. Let 

n = -I- nib^ 


represent the base-6 expansion of n and denote di := di{n) = ni to be the i-th digit of n in base b. We shall 
say that m is digitally less than or equal to n if rui < ni for all j = 0, 1, ..., A — 1. In that case, we shall 
write m :< n. We are now ready to state our result. 
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Theorem 2. Let n be a non-negative integer. Then 

iV-l 


n 

i=0 


x + y + di(n) - 1 

di{n) 


= E 

0<m-<n 


^ ^ 'x + di{m) — 1\ /y + di{n — m) — 1 


n 


di{m) 


n 


di{n — m) 


(5) 


Let yij{n) := {n) denote the multiplicity of the digit j > 0 in the base-6 expansion of n, i.e., 

yj{n) = |{i : d,{n) = j}\. 

As a corollary, we obtain 

Corollary 3. Let n be a non-negative integer. Then 

b-l 




3=0 


J 


0<m-<n 


n 

i=i 


X J — -p-r (y J — 


J 


n 


J 


Observe that if 6 = 2, then Corollary [3] reduces to Theorem [T] 

The source behind Theorem [1] is a one-parameter subgroup of Sierpinski matrices, investigated by Callan 
in [3], which encodes the digital binomial theorem. To illustrate this, define a sequence of lower-triangular 
matrix functions Sn (x) of dimension 2^ x 2^ recursively by 


Si{x) = [ ^ p ): Sn+i{x) = Si{x) ® SNix), 


( 6 ) 


where 0 denotes the Kronecker product of two matrices. For example, S 2 {x) and 53 ( 0 ;) can be computed as 
follows: 


S 2 {x) = Si{x) 0 Si(x) = 




/ 

1 

0 

0 

0 ^ 





1-Si{x) 

0- Si{x) \ _ 


X 

1 

0 

0 





X ■ Si{x) 

1- Si {x) J ~ 


X 

0 

1 

0 







V 

x^ 

X 

X 

1 







/ 

1 

0 

0 

0 

0 

0 

0 

0 




X 

1 

0 

0 

0 

0 

0 

0 




X 

0 

1 

0 

0 

0 

0 

0 

1 • S2{X) 

0 • 52(2:) \ _ 


x^ 

X 

X 

1 

0 

0 

0 

0 

X ■ S2{x) 

1 • S2ix) J ~ 


X 

0 

0 

0 

1 

0 

0 

0 




x^ 

X 

0 

0 

X 

1 

0 

0 




x^ 

0 

X 

0 

X 

0 

1 

0 



1 

x^ 

x^ 

x^ 

X 

x^ 

X 

X 

1 


Sslx) = Si(x) 0 S 2 (x) = 


Observe that when a: = 1, the infinite matrix S = lim^v-^oo 'Siv(l) is of course Sierpinski’s triangle. 

In |3], a formula for the entries of Sn{x) = {aN{j,k,x)), 0 < j, fc < 2^ — 1, is given in terms of the 
sum-of-digits function: 


unU, k,x) = 


k) ^ if 0 < fc < j < 2^ — 1 and (fc, j — k) are carry-free 


0 , 


otherwise. 


(7) 


Moreover, it was proven that Sn{x) forms a one-parameter subgroup of SL{2^, 
real matrices with determinant one. Namely, we have 

SN{x)SN{y) = Sn{x -L y). 

If we denote the entries of SN{x)SN{y) by tiq{j, k), then the equality 

2~-l 


i.e., the group of 2^ x 2^ 
( 8 ) 


iN{jjk)= ^ aNijji,x)aN{i,k,y) = aN{jjk,x y) 


2=0 
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corresponds precisely to Theorem [T] with j = 0 and k = 0. For example, if = 2, then ([ 5 ]) becomes 


/1 

0 

0 

0 ^ 


/1 

0 

0 



/ 1 

0 

0 

0 \ 

X 

I 

0 

0 


y 

1 

0 

0 


x + y 

1 

0 

0 

X 

0 

1 

0 


y 

0 

1 

0 


x + y 

0 

1 

0 

V 

X 

X 

1 ^ 


\ y^ 

y 

y 

1 ^ 


\ (x + y)^ 

x + y 

x + y 

1 / 


The rest of this paper is devoted to generalizing Callan’s construction of Sierpinski matrices to arbitrary 
bases and considering two applications of them. In Section 2 , we use these generalized Sierpinski matrices 
to prove Theorem [21 In Section 3 , we demonstrate how these matrices arise in the construction of Prouhet- 
Thue-Morse polynomials defined in m- In Section 4 , we describe a group presentation in terms of generators 
defined through these matrices and show that these generators satisfy a relation that generalizes E. Ferrand’s 
result in [3]. 


2. Sierpinski Triangles 


To prove Theorem |21 we consider the following generalization of the Sierpinski matrix Sn{x) in terms of 
binomial coefficients. Define lower-triangular matrices Sb^N{x) of dimension x recursively by 



f 

1 

0 

0 

... o\ 





(?) 

1 

0 

... 0 



Sb,lix) = 


(?) 

1 

... 0 

II 

+ 

1 1 

1 

if0<fc<j <6—1 

[ 


/a:+b—3'\ 

V b-2 J 

fx+b—4\ 

V 6-3 / 

... l) 

1 0, 

otherwise. 



and for > I, 

Sb,N+l{x) = Sb,lix) (g) Sb,Nix). 

Example 4. To illustrate our generalized Sierpinksi matrices, we calculate Ss^ii^x) and S3^2{x): 





1 0 

0 \ 








S3,l(x) 

= 


(2 1 
:t^) (?) 










[ 

J 








S3,2(x) 

= S3,l 

(x) o S'3,1 

(x) 










f 

1 

0 

0 

0 

0 

0 

0 

0 




(?) 

1 

0 

0 

0 

0 

0 

0 

0 




(X) 

(?) 

1 

0 

0 

0 

0 

0 

0 




0 

0 

1 

0 

0 

0 

0 

0 


= 


(?)(?) 

(?) 

0 

(?) 

1 

0 

0 

0 

0 




(?)(%^^) 

(7^) 

mil) 

mm 

(?)(?) 

(?) 

(^r) 

(?) 

1 

0 

0 

0 




0 

0 

(?) 

0 

0 

1 

0 

0 




(%^^) 

mil) 

0 

(?)(?) 

(?) 

0 

(?) 

1 

0 



[ 

(^r) 

(?)(“r) 

(?)(?) 

(?) 

(^r) 

(?) 

1 / 


We now generalize Callan’s result for Sjsf{x) by presenting a formula for the entries of St^wix) (see [ 7 ] for 
a similar generalization of Callan’s result but along a different vein). 

Theorems. Let aAr(j, fc) := ajv(j, fc, x) denote the {j,k)-entry of Sb,N(x). Then 

\ ri"Arir‘)- cSir’)’ 

ajv(jA)=< ( 9 ) 

[ 0, otherwise. 

where j — k = dob^ -f dib^ -f ... -f dLb^ is the base-b expansion of j — k, assuming j > k. 
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Proof. We argue by induction on N. It is clear that holds for Sb,i{x). Next, assume that m holds 
for Sb,Nix) and let aN+i{j,k) be an arbitrary entry of Sb,N+iix), where pb^ < j < {p + l)b^ — 1 and 
qb^ < k < {q + l)b^ — 1 for some non-negative integers p,q £ { 0 , 1 ,..., 6 — 1}. Set j' = j — pb^ and 
k' = k — qb^. We consider two cases: 

Case p < q. Then j <k and a7v+i(j, fc) = 0 • k') = 0 . 

Case 2. p> q. Then j > k and 

aN+i{j,k) = ^ ^ ^'\aN{j',k'). ( 10 ) 

V P-9 J 


Let j — k = dob° + dib^ dNb^, where djv = P — 9- Then j' — k' = -I- dib^ dn-ib^ By 

assumption, 


aN{j',k') = 


(.+*-.) (■+*-!),.. if 0 < r < / < - 1 and t' /, 


( 11 ) 


0 otherwise. 

Since k ^ j ii and only if k' ^ j', it follows from (HHD and m that 

CT""') (“h'r‘) ^ ■ ■ ('dr‘) if 0 < * < i < <>~+‘ -1 and t d j. 

0 otherwise. 


aN+iij,k) = 


( 12 ) 


Thus, (O holds for Sb,N+i- 


□ 


Next, we show that Sb,N{x) forms a one-parameter subgroup of SL{b^,'R). To prove this, we shall need 
the following two lemmas; the first is due to Gould [B] and the second follows easily from the first through 
an appropriate change of variables. 

Lemma 6 (Gould [B]). 

^ /'x + k\/'y + n — k\ /x + y + n + l 
k )\ n-k ) ~\ n 

Proof. In [B] , Gould derives (ITBl) as a special case of a generalization of Vandemonde’s convolution formula. 
We shall proof ( 1131 ) more directly by presenting two different proofs, one using a combinatorial argument 
and the other using the beta function. 

Combinatorial argument. Let A, B, and C = { 0 , 1 ,..., n} denote three sets containing x, y, and n -I- 1 
elements (all distinct), respectively, where n is a non-negative integer. For any non-negative integer fc, define 
Ak = AU{ 0 ,..., fc — 1 } and Bk = B\j{k+ 1 ,..., n}. Then given any n-element subset S of AUi?UC, there 
exists a unique integer fcg in C — S, called the index of S with respect to A and B, such that [S' n Ak^ |= ks 
and \S l~l Bfcg \= n — ks. To see this, define SA = Ar\S,SB = BnS, and T = C — S. We begin by deleting 
IS'aI consecutive elements from T, in increasing order and beginning with its smallest element, to obtain a 
subset T' . We then delete IIS's! consecutive elements from T', in decreasing order and beginning with its 
largest element, to obtain a subset T", which must now contain a single element denoted by ks. It is now 
clear that jS” D ^fcs 1 = and jS” n Bks\= n — ks- 

To prove (fTBl) . we count the n-element subsets S' of A U i? U C in two different ways. On the one hand, 
since U B U C\= x + y + n+ 1 , the number of such subsets is given by Qn the other hand, we 

partition all such n-element subsets into equivalence classes according to each subset’s index value. Since 
|S n ^fc|= k and \S r\ Bk\= n — k, it follows that the number of n-element subsets S having the same index 
k is given by and total number of n-element subsets is given by 

fx + k\fy + n - k^^ 

^ I k j I n-k )' 

fe=0 ^ z ^ z 

Lastly, we equate the two answers to obtain m- 
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Analytic argument. Recall that the beta function is defined as B{x,y) = r(a:)r(j/)/r(a; + y), where r(a;) is 
the gamma function. If x is a non-negative integer, then r(x -I- 1 ) = x! and 


B{x,y) 


(x- l)!(y- 1)! 
{x + y- 1)! 


For convenience, we shall write x! to represent r(x -I- 1 ) even when x is not non-negative. We now divide 
m hy n\x\y\ /{x + y -\- n + 1 )\ to obtain the following identity: 



B[x + k + l,y n 


fc -I- 1 ) = B{x -I- 1 ,?/ -I- 1 ). 


But (full is a generalization of the following classical property of the beta function: 

R(x, y -l) + B{x -l,y) = B{x -l,y- 1 ). 

It is now straightforward to prove ( 1141 ) using an induction argument. 


( 14 ) 


( 15 ) 

□ 


Lemma 7 . Let p and q be positive integers with q < p. Then 


E 

v—q 


x+p —V —l\fy + v —q — 1 \ fx + y+p — q —1 


p — V 


v — q 


p-q 


Proof. Set w = V — q. Then dni) can be rewritten as 


E 

w—O 

which follows from Lemma [51 
Theorem 8. For all N gN, 


x+p —q —w —l\fy + w — 1 \ fx + y+p — q —1 


p — q — w 


p-q 


Sb,Nix)Sb,Niy) = Sb.N^x + y). 


( 16 ) 


□ 

( 17 ) 


Proof. We argue by induction on N. Lemma [6] proves that (ITTI) holds for Sb.i{x). Next, assume that (ED 
holds for Sb,Nix). Let tN+iij,k) denote the (j, fc)-entry of Sb,N+i{,x)Sb,N+i{y)- Then 

b^+^-i 

tN+iij,k)= ^ aN+iij,m,x)aN+i{m,k,y) 
m—0 

b-l &"-l 

= ^ E +'>’,x)aN+i{vb’^ +r,k,y). 

v=0 r—0 

As before, assume pb^ < j < {p + l)b^ — 1 and qb^ < k < {q + l)b^ — 1 for some non-negative integers 
p,q € { 0 , 1, ..., 6 — 1}. Set j' = j — pb^ and k' = k — qb^. Again, we consider two cases: 


Case 1 : j < k. Then ajv+i(j, fc, x + y) = 0 by definition. On the other hand, we have 

j 

tN+iij,k) = '^ aN+i{j,rn,x)aN+i{m,k,y) + ^ ajv+i(j, m, x)ajv+i(m, fc, y) 

m—0 rn—j-\-l 


j 

= E “JV+i(j, W,x) • 0 + ^ 0 • aAr+i(m, fc, y) 

m—0 m=j-\-l 

= 0 


and thus ED holds. 


Case 2 : j > fc. Since Sb,N+i{x) = Sb,i{x) <S> Sb.Nix), we have 


if j > vb^ + r, 
if j < vb^ + r. 
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Similarly, we have 


uN , „ _ f ^)aN{r,k',y) ii k < vb^ + r, 


aN+iivb^ +r, fc, y) = 


if A: > vb‘^ + r. 


It follows that 


tN+i{j,k) = J 2 Y 1 

v—q r—0 


x-\-p — V — l\jy-\-v — q —1 
p — V J \ v — q 


X p — V — l\ fy -\- V — q — 1 
p — V J \ v — q 


aN{j',r,x)aN{r,k',y) 
^ aNU\r,x)aN{r,k',y) 


aN{j',k',x + y) 


_v—q N ^ /V ^ ' J r—0 

''x + y + p — q — V 
p-q 
= aN+i{j,k,x + y), 

where we have made use of the inductive assumption and Lemma 0 This proves that ( 1171 ) holds. 

As a corollary, we obtain Theorem [ 2 J which we now prove. 

Proof of Theorem\^ Let j = n and k = 0 . Then the identity 

b^-i 


□ 


aN{j,m,x)aN{rn,k,y) = aN{j,k,x + y), 


m—0 


which follows from (HU, is equivalent to dS]). 


□ 


We end this section by describing the infinitesimal generator of Sb,N{x). Define Xb.i(x) = {xj,k) to be a 
strictly lower-triangular matrix whose entries Xj,k are given by 

x/{j-k), iij>k + l 
0 , otherwise. 


Xj,k = 


( 18 ) 


For > 1 , we define matrices 

Xb,N+lix) = Xb^i{x) © Xb^N{x) = Xb^i{x) 0 IbN + Ib0 Xb^NiV^ 

where © denotes the Kronecker sum and Ib^ denotes the b^ x b^ identity matrix. Observe that Xb,i{x) has 
the following matrix form: 


Xb,i{x) = 


X 

0 

0 

... u 1 

... 0 

xl2 

X 

0 

... 0 

xfi 

x/2 

X 

... 0 

\ a:/(&- 1) 

xj {b — 2) 

xj {b — 3 ) 

... o) 


( 19 ) 


The following lemmas will be needed. The first states a useful identity involving the unsigned Stirling 
numbers of the first kind, c(n, k), defined by the generating function 

n 

x{x + l)...(a: + n — 1 ) = c(n, k)x^ . 

k=0 

It is well known that c(n, k) counts the number of n-element permutations consisting of k cycles. 

Lemma 9. Let I and n be positive integers with I > n. Then 


Z —n+1 


" /l\ 

(j — 1)! ( . jc(l — J, n — 1) = nc(Z, n). 


( 20 ) 
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Proof. We give a combinatorial argument. Let A = {1, 2, 1}. We count in two different ways the number 
of permutations tt = cria 2 ■ • ■ (Jn of A consisting of n cycles where we distinguish one of the cycles aj of 
TT. On the one hand, since there are c{l,n) such permutations tt and n ways to distinguish a cycle of tt, it 
follows that the answer is given by nc{l,n). On the other hand, we can construct tt by first choosing our 
distinguished cycle ui consisting of i elements. The number of possibilities for a\ is (*)(* — 1)! since there 
are (*) ways to choose i elements from A and (i — 1)! ways to construct a cycle from these i elements. It 
remains to construct the remaining cycles (T 2 , • • •, cr„, which we view as a permutation tt' = ■■■ <Jn on I — i 

elements consisting of n — 1 cycles. Since there are c{l — i,n — 1) such possibilities for tt', it follows that the 
number of permutations tt with a distinguished cycle is given by 

l — n + l y,\ 

Equating the two answers yields as desired. □ 


Lemma 10. Let n be a positive integer with 1 < n < b — 1. Then 


= {XnU,k)), 


( 21 ) 


where the entries Xnij, k) are given by 


Xn{j,k) 


n! 

O-fc)! 


c(j — fc, n)a:", 


0 , 


if j > k + n 
otherwise. 


( 22 ) 


Proof. We argue by induction on n. It is clear that (l22l) holds when n = 1. Suppose n > 1. li j < k + n, 
then Xn{j, k) = 0 because is a power of strictly lower-triangular matrices. Therefore, assume j > k + n. 
Then 


h—1 j—n+1 

Xn{j,k) = '^Xn-i{j,i)xi{hk) = ^ Xn-i{j,i)xi{hk) 

2—0 2 — fc +1 

j-n+1 / . N. . 




(j -*)! 


i — k 


(n — 1)! a;’ 

Tl 




X (*“!)■ .]c{l-i,n-l), 


where I = j — k. It follows from Lemma [5] that 


Xn{j,k) = 


(n — 1)! X 


n\ 


■nc{l,n)= J - k,n)x^‘ 


as desired. 


□ 


Lemma 11. We have 


exp{Xb^i{x)) = Sb,i{x). 


(23) 
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Proof. Denote the entries of exp(Xb,i(x)) by ^(j, k). It is clear that ^{j, k) = 0 for j < k and ^{j, k) = 1 for 
j = k since is strictly lower triangular. Therefore, assume j >k + l. Since = 0 for n > 5, we have 


Thus, (1^ holds. 




Xn(j, k) 
^ nl 


1 


b-l 


{j-ky. 


^ c{j - k, n)a 


n—1 

x{x + 1) ■ ■ ■ {x + j — k — 1) 


(j - ky 

fx + j- k-l^. 

I J-k J 

ai{j,k). 


□ 


Theorem 12. Let N he a positive integer. Then 

exp{Xb^Nix)) = Sb,N{x). (24) 

Proof. We argue by induction on N. Lemma (ED shows that ((Ml) is true for N = 1. Then since exp(A©i3) = 
exp(y4) © exp(i3) for any two matrices A and B, it follows that 

exp(Xb,Ar(a;)) = exp(X(,,i(x) © Xb,N-i{x)) = Sb,lix) © Sb,N-iix) = Sb,N{x), 

which proves (IMl) . □ 


3. Prouhet-Thue-Morse Polynomials 


In this section we demonstrate how the generalized Sierpinski matrices Sb.Ni^) arise in the study of 
Prouhet-Thue-Morse polynomials, first investigated by the author in m- These polynomials were used in 
the same paper to give a new proof of the well-known Prouhet-Tarry-Escott problem, which seeks b > 2 sets 
of non-negative integers So, Si, ..., Sb-i that have equal sums of like powers up to degree M > 1, i.e., 

^ ^ n"* = • • • = 

n£So n€Si n€Sb-i 

for all m = 0,1,..., M. In 1851, E. Prouhet El gave a solution (but did not publish a proof; see Lehmer [8]) 
by partitioning the first 6 ^+^ non-negative integers into the sets So,Si,..., Sb-i according to the assignment 


^ ^ ^Ub{n) • 

Here, Ub{n) is the generalized Prouhet-Thue-Morse sequence, defined as the residue of the sum of digits of 
n (base b): 

d 

Ub(n) = nj mod b, 
i=o 

where n = ndb‘^ + -h nob° is the base -6 expansion of n. When 6 = 2, u{n) := U2{n) generates the classical 

Prouhet-Thue-Morse sequence: 0,1,1, 0,1, 0,0,1, _ 

Let A = (ag, ai,..., Ob-i) be a zero-sum vector, i.e., an ordered collection of 6 arbitrary complex values 
that sum to zero: 


Oo T T ■ ■ ■ T Ob—1 — 0- 

We define Fn{x;A) to be the Prouhet-Thue-Morse (PTM) polynomial of degree 6 ^ — 1 whose coefficients 
belong to A and repeat according to Ub{n), i.e.. 


Fn[x-,A) = ^ 

n—0 
8 




(25) 





In the case where 6 = 2, ao = 1, and oi = —1, we obtain the classic product generating function formula 

AT 2 "+i-l 

n(i-^^'")= E (-1)“^"^^”- (26) 

m—O n—0 

Lehmer generalized this formula to the case where A= consists of all 6 -th roots of unity 

with uj = The following theorem, proven in m, extends this factorization to F{q{x]A) for arbitrary 

zero-sum vectors. 


Theorem 13 1[10]'). L^t. &6 d pOSli/tUG ’Ifltc.QG.T CLTid Q, ZGTO SUTTl VGCtOT, ThGTG GXZSts (X polyThOTThtCll 
such that 

N-l 

Fn{x;A) = Pn{x)Y[{1-x^’"). (27) 

m—0 


Theorem [T3] is useful in that it allows us to establish that the polynomial Fiy{x, A) has a zero of order N 
at a; = 1, from which Prouhet’s solution follows easily by setting N = M + 1 and differentiating Fpf{x; A) m 
times (see m)- 

We now derive formulas for the coefficients of PnIx) in terms of generalized Sierpinski triangles. Towards 
this end, let 

Pn{x-,Cn) = E 

n—0 

denote a polynomial whose coefficients are given by the column vector cat = (cq, ..., ChN _ i )'^ . Also, let 

(^Ub{0); ^Ub(l) ; * ■ ■ ; “1) ) 


be a column vector consisting of elements of A generated by the PTM sequence Ub{n). Next, define a 
sequence of x matrices Mjv recursively by 


and for > 1 , 


Mn+1 


Ml = 


Ml ®Mn = 


/ 1 0 

0 .. 

0 




-1 1 

0 .. 

0 

0 



0 0 

0 .. 

. -1 

1 ) 



( Mn 

Oat 

Ojv 


On 

Oat \ 

On 

-Mn 

Mn 


On 

On 

\ Oat 

On 

On 


—Mn 

Mn ) 


(28) 


where Mi ® M^ denotes the Kronecker product. The following theorem establishes a matrix relationship 
between the vectors and cat. 


Theorem 14. Let A = (oq, ..., a&-i) be a zero-sum vector. The polynomial equation 

N-l 

Fn{x-A) = Pn{x-Cn) n (1 -^'''") ( 29 ) 

m—0 

is equivalent to the matrix equation 

aAT = MjvCjv (30) 

together with the condition c„ = 0 for any n that contains the digit b — 1 in its base-b expansion, where 
0 <n <b^ — 1. 


To prove Theorem lldl we shall need the following lemmas, which we state without proof since their results 
are easy to verify. 
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Lemma 15. Let A = (ao, ai, ..., ab-i) and C = {cq, ■ - ■, Cb-i}- Then the polynomial equation 


6-1 

n—0 

is equivalent to the system of equations 



(l-x) 


ao = Co 

ai = —Co + Cl 


flb-l = —Cb-2 + Cb_l 


together with the condition Cb-i = 0 . 

Lemma 16. The system of equations in Lemma [751 can he expressed in matrix form as 


ai = MiCi- 


Proof of Theorem We argue by induction on N. It is clear that (1501) holds for TV = 1 since the polynomial 
equation Fi{x; A) = Pi(x; C'i)(l —x) is equivalent to ai = MiCi because of Lemmas fTSl and fTBl Next, assume 
that (l30l) holds for case N. We shall prove that (l30l) holds for case + 1. Define 

b^-i 

Pn{x]Cn{p)) = 

n—0 

to be a polynomial with coefficient set Cn{p) = {cpb-w,..., C(p+i)f,N-i}- We then expand the right-hand side 
of (l?7l) for case -b 1 as follows: 


N 

Pn+i{x\Cn+i) (1 

m=0 


) = [-PAr(^; CivlO)) + ... + Pn{x'-, CNib — 1))] 


'N-1 


_m—0 


(1 


— [QNix; C'jv(O)) -b ... -b Qn{x', Cj\[{b — 1))](1 — x^ ) 

= QNix; C7v(0)) -b [QNix; C'Ar(l)) — x^ QNix; C'iv(O))] -b ■ • ■ 

lN iN 

+ [QNix;CNib-l)) - x^ QNix;CNib-2))]-x^ (3jv(x;C'A r(& - 1)), 


,])j 

where we define QNix; Cjsiip)) = PNix; C' 7 v(p))(l — x° ). Equating this result with 

bJV + l_i 

FN+iix;A)= a„^(„)x" = F/v(x; Ajv(O))-b ...-b E/v(x; — 1)), 

n—0 


where 


FNiXf Ajsbipy) ^ ^ ^u{n+pb^)^ ^ : 

n—0 


leads to the system of polynomial equations 


F/v(x; Ajv(O)) = QNix; C'jv(O)) 

lN 

FNix; AAr(l)) = QNix; C'Ar(l)) — x^ QNix; C'7v(0)) 


FNix; ANib — 1)) = QNix; C'jv(& — 1)) — x^ QNix; CNib — 2)) 
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together with the condition Qn{x; CN{b — 1)) = 0, i.e. c„ = 0 for all {b — l)b^ < n < 6 ^+^ — 1. It follows 
from Lemmas [15] and [T 6 | that this system is equivalent to the system of matrix equations 

aAr(O) = MnCn{0) 

aAr(l) = —MnCn{0) + MnCn{1) 


aN{b — 1 ) = —MjvCjv(& - 2 ) + M^CNlb — 1 ), 

where the first matrix equation by assumption satisfies the condition c„ = 0 for any 0 < n < b^ — 1 
that contains the digit 5 — 1 in its base -6 expansion. This in turn is equivalent to the matrix equation 
a^r+i = Mjv+iCjv+i together with the condition that c„ = 0 for any 0 < n < 6 ^+^ — 1 that contains the 
digit 6 — 1 in its base -6 expansion. This establishes the theorem for case -|- 1 and completes the proof. □ 


Lemma 17. The matrix Mjq has inverse Sm 

= 

Mm\ 

where Sm is given recursively by 



/ 1 0 

0 .. 

■ ^ ) 






I 1 

0 .. 

. 0 




Si = 







(31) 


U 1 

1 .. 

• iJ 




and for N > 1, 


/ Sm 

Oat 

On 

... Oat \ 





Sm 

Sm 

Om 

... Oat 



Sn+1 = Si iSi Sm -- 


\ Sm 

Sm 

Sm 

... Sm j 


(32) 

Thus, if aN = MmCm, then 









Cm = Sm^n- 




(33) 


Proof. It is straightforward to verify directly that Since S'jv+i = S'i( 8 >S'jv and Mjv+i = Mi®Mjsi, 

it follows that Sn = □ 


Observe that Sm = Sb^Ni^) is the generalized Sierpinski triangle defined in the previous section. We now 
present a formula for the coefficients c„ in terms of the elements of A. 

Theorem 18. Let A = {ag,..., ab-i} be a zero-sum vector. Then the polynomial equation P29\) has solution 
Pn{x', Cat) whose coefficients Cn, 0 < n < b^ — 1, are given by 

^ ^ ^Ub{k)i 
fce/|,(n) 

where Ib{n) = {fc S N : fc ^ n}. Moreover, c„ = 0 for all n whose base-b expansion contains the digit 6—1. 

Proof. From Theorem (5] we know that the non-zero entries in the n-th row of Sn, which are all equal to 1, 
are located at (n, k) where k < n. Formula (IMl) now follows from (1551) . It remains to show that (1551) yields 
Cn = 0 for all n whose base -6 expansion contains the digit 6—1. Let n = ngb^ -(-... -i-n^b^ n^-i bN-l 

where = b — 1. Denote Ib{n-,L) to be the subset of Ib{n) consisting of integers whose base -6 expansion 
has digit 0 at position L, i.e., 

Ib{n-, L) = {k £ N : k n, k = kob^ k^b^, kL = 0}. 

Then using the fact that >1 is a zero-sum vector, i.e., og -I- ... -I- Ob-i — 0, we have 

Cn = ^ aui,(fc) + ^ C,n^{k+bV) + ■ ■ ■ + ^ ani,(k-b{b-l)bV)) 

fc€/b(n;L) kGlb{n\L) k^Ib{n\L) 

= X! [®«6(fe) + ®(«fc(fe)-|-l)6 + ■ ■ ■ + ®(«fc(fc)-|-(b-l))6] 
keiIb(n;L) 

= 0 

□ 


as desired. 
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Next, we specialize Theorem [T51 to base 6 = 3. Define w{n) to be the sum of the digits of n in its base-3 
representation modulo 2 . 

Corollary 19. Suppose 6 = 3. Let n gN be such that its base-3 representation does not contain the digit 2. 
Then 

c„ = (-l)-Wa„ 3 ( 2 „). (35) 

Proof. Let n = be the base-3 representation of n with no digit equal to 2 and leading 

digit = 1. We argue by induction on N. It is clear that (1551) holds when N = 1 since from ((Ml) we 

have 

Cl Uq ^2 ( 1) ^ ^ ^U3{2n) • 

Next, assume dSSl) holds for a given N. To prove that (1551) holds for iV -|- 1, we consider two cases by 
decomposing n = n' -|- 3^. 


Case 1: = 0 for alH = 0,..., — 1. Then n = 3" and 

^ ^ ^u^{k) ^113(0) “f ^u^{n) ^0 “f ^2 

kel3{n) 

_ f 1 \w(n) 

( 1 ) 0'ii3(2n) • 

Case 2: n/, = 1 for some 0 < L < — 1. Then set n" = n' — 3^. It follows that 

Cn — ^ ' ®ii3(fe) ^ ' ®U3(fe) T ^ ' ®ii3(fc-|-3''^) 

fe£/3(n) fcg/3(n,L) kGl3{n,L) 

feG/3(ra') kelsin') 

— 'y ' Ou3(fe) + ^ ^ ®ii3(fe-|-3^) T ^ ^ Oii3(fc-|-2-3^) T ^ ' ^U3{k) ~ 'y ' 0,U3{k) 

kalsin') kGl3{n' ,L) kGl3{n' ,L) kel3{n' ,L) kC:l3(n' ,L) 

~ ^ ' fl«3(fe) T ^ ^ [®«3(fc) ®(«3(fc)-|-l)3 ®(M3(fc)-|-2)3] ~ ^ ^ fl«3(fc) 

feG/3(n') kelsin' ,L) kelsin") 

( Ij ^it3(2n^) ( Ij ^it3(2n") 

= (-1)™^” ^[a„3(2„') -I- a(„3(2„')_2)3] = ^[““^(«3(2n')-l)3] 

— ( —Ij <^{it3(2n0+2)3)) — V“l) ^ii3(2n'+2-3^) 

— ( 1 j ^U3(2n) • 

Thus, ((551) holds for A^ -I- 1. 


□ 


4. Group Generators and Relations 

In this section, we describe group generators and relations defined by the matrices Sjy and Mjy. Recall 
that Sjv and Mjy are matrices of dimension 6 ^ x 6 ^ for a given base 6 . Define Tjv = to be the transpose 
of M]y and Un = SnTjsi, Vn = T^Sn. The following lemma gives a recursive construction of Cat and Vat. 

Lemma 20. We have 

Un+ 1 = Ui® Un 
Vn+1 = Vi <8) V/v 

Proof. The result follows from the mixed-product property of the Kronecker product. We demonstrate this 
for Dat+i: 

Un+1 = Sn+iTn+i = {Si ® Sn){Ti ® Tn) = (SiTi) ® {SnTn) = Di ® Un- 

□ 


The calculation is the same for Vw+i. 
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Observe that Un = {uij) and Vn = (vij) are skew-triangular and that Vn is the skew-transpose of Un, 
i.e., Vij = Ub-i-j,b-i-i ior i,j = 0,1,... ,b — 1. Here are some examples of Un and Vn when b — 2: 









(1 

-1 

-1 

1 

-1 

1 

1 

-1 \ 








1 

0 

-1 

0 

-1 

0 

1 

0 



( 1 

-1 

-I 

1 ^ 


1 

-1 

0 

0 

-1 

1 

0 

0 

^ -1 \ 

, U2 = 

1 

0 

-I 

0 

, Uz = 

1 

0 

0 

0 

-1 

0 

0 

0 

V 1 0 ) 

1 

-1 

0 

0 

1 

-I 

-1 

1 

0 

0 

0 

0 



1 1 

0 

0 

0^ 


1 

0 

-1 

0 

0 

0 

0 

0 






1 

-I 

0 

0 

0 

0 

0 

0 








11 

0 

0 

0 

0 

0 

0 

0 / 









/ 

0 

0 

0 

0 

0 

0 

0 

\ 









0 

0 

0 

0 

0 

0 

1 

1 



/ 0 

0 

0 

1 ^ 



0 

0 

0 

0 

0 

1 

0 

1 

/ 0 -1 ^ 

, 1^2 = 

0 

0 

-I 

-1 

, 1^3 = 


0 

0 

0 

0 

-1 

-1 

-1 

-I 

(l 1. 

0 

-1 

0 

-1 


0 

0 

0 

1 

0 

0 

0 

1 



11 

1 

1 

1 ) 



0 

0 

-1 

-1 

0 

0 

-1 

-I 







0 

-1 

0 

-1 

0 

-1 

0 

-I 








V 

1 

1 

1 

1 

1 

1 

1 

1 / 


The next lemma establishes that the eigenvalues oi Un and Vn are (6 -I- l)-th roots of 1 or — 1. 


• (37) 


Lemma 21. The set of eigenvalues of Ui and Vi are exactly the same and consist of all roots of the 
polynomial equation 

-1 + r -r"^ + ... + = 0. 

Proof. Let r be a root of —1 -I- r — -I- ... -I- (—= 0. We claim that r is an eigenvalue of Ui with 

eigenvector v = (ui,..., Vb)'^, where 

0 , k = b. 

Denote w = {Ui — rlb)v = {wi,..., Wb}. It suffices to show w = 0. Assume b = 2. We have 

= (1 — r)r — l = —l-\-r —r^ = 0 
W 2 = r — r = 0. 

For b > 2, we have 

wi = (1 — r)r — {r — r^) = 0 

W 2 = r - r{r - r^) - {r - r'^ + r^) = 0 



Wb-i = r — r(r — r'^ + r^ + ... + r’^ ^) — 1 = —1 -I- r — -I- ... -I- (—= 0 

Wb = r — r = 0. 

Thus, w = 0. It can be shown by a similar argument that the eigenvalues of Vi are exactly the same as 
those of C/i. □ 


Theorem 22. The matrices Un and Vn satisfy the relation 

Un^" = Vn'"" = 

where Un is the b^ x b^ identity matrix. 


(38) 


Proof. We shall only prove (IMll for Un since the proof for Vn is the same. We argue by induction on N. 
When = 1, we know from Lemma[2T]that all eigenvalues of Ui are roots of —1-l-r —r^-|-.. .-|-(—= 0. 
It follows that every eigenvalue r satisfies = (—1)^+^ and since they are all distinct, the corresponding 
eigenvectors are all linearly independent. Thus, U^'^^ = (—I)^+^/b. 
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Next, assume that (1551) holds for Un-i- It follows from the mixed-product property of the Kronecker 
product that 

= (C/i 0 ® 

Thus, (1551) holds for □ 

We note that for h = 2, Ferrand proved in [?] that the matrices Sn and Tm satisfy the 3-strand braid 
relation 

SnTnSn = TnSnTm- (39) 

We give an alternate proof of (1551) based on Theorem 1551 Define Qn = SnTnSn and Rn = TmSnTn- Then 
QnRn = Uff = (—= (—l)'^J 2 iv because of (1551) . Moreover, we claim that = Rff = {—1)^12^. 
This follows by induction on N, which we demonstrate for by again using the mixed-product property 
of the Kronecker product: 

Q'n = (SnTnSIiTnSn) 

= (SiTiSfTiSi) (g) iSN-iTN-iSlf_^TN-iSN-i) 

= (-l)^J2iV 

Thus, QnRn = which implies the braid relation Qjv = Rn- However, we find that the braid relation 
fails to hold for & > 2. 
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